Dynamics and decoherence in the central spin model using exact methods 
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The dynamics and decoherence of an electronic spin-1/2 qubit coupled to a bath of nuclear spins 
via hyperfine interactions in a quantum dot is studied. We show how exact results from the integrable 
solution can be used to understand the dynamic behavior of the qubit. It is possible to predict the 
main frequency contributions and their broadening for relatively general initial states analytically, 
leading to an estimate of the corresponding decay times, which are related to T\ of the electron. 
Furthermore, for a small bath polarization, a new low-frequency time scale is observed. 
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One of the most challenging tasks in both theoretical 
and experimental studies of quantum information pro- 
cessing is the understanding and control of decoherence 
due to the interaction with the environment. A proto- 
typical example for the loss of quantum information is 
an electron trapped in a quantum dot, surrounded by nu- 
clear spins. On short time scales up to 1ms, ^ the Heisen- 
berg exchange resulting from the hyperfine interaction 
between the electron and the nuclei dominates, before 
spin-orbit coupling or dipole-dipole-interactions between 
the bath spins become effective. ^''^ Indeed, impressive ex- 
perimental progress has been made over the recent years 
to observe and control oscillations of single electron spins 
coupled to bath spins by various techniques. ^^^ In a typi- 
cal magnetic resonance experiment, for example, one ap- 
plies a static magnetic field creating a Zeeman splitting of 
the spin's energy levels, plus a resonant or near resonant 
alternating field transverse to the Zeeman field. Spin ro- 
tation angles can then be controlled e.g. by the duration 
of the transverse field pulse. The decoherence processes 
on which we focus here, however, are not directly linked 
to the applied fields, but are instead caused mainly by 
the Heisenberg coupling to the nuclear bath spins, which 
is given by Hamiltonian 
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Nt 

E 



A J So ■ Sj , 



(1) 



where Nf, is the number of bath spins and So is the central 
electron spin. This so called central spin model is one of 
the most studied theoretical models for decoherence, ^"^^ 
which can also be treated by exact methods. ^^"^^ In this 
paper we are interested in an analytic understanding 
of the detailed decoherence- induced dynamics of {So){t) 
due to the coupling to the spin bath. We therefore do not 
consider any external magnetic fields, but instead spec- 
ify the initial overall polarization of the system, which is 
conserved. The central spin Sq is initially assumed to be 
in its down state, independent of the bath spins. This 
initial product state gets entangled by the exchange in- 



teractions, leading to the decoherence of the central spin. 
The couplings Aj > in Eq. (1) are proportional to 
the square of the electronic wave function at the position 
of the nucleus j. The methods we apply in the follow- 
ing do not depend on the special choice of the Aj . For 
definiteness, we assume a Gaussian distribution with the 
site index (distance) ji'^'^i 
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which allows an easy control over the two relevant char- 
acteristics of the distribution of Aj, namely the mean 
value xi = j^X^i^j ^^'^ ^^^ degree of inhomogeneity 
as parametrized by B. Generally the results are largely 
insensitive to the overall shape of the distribution, (e.g. 
in case a higher dimensional site index is used) as long 
as the mean xi and the degree of inhomogeneity are the 
same. 

For homogeneous couplings, Aj = A\fj, a non-trivial 

time scale T'^ A~^Nf^ has been identified using exact 
methods, ^^'^'^ which can be interpreted as a decoherence 
time. A number of authors'^'-'^^"^* have studied the in- 
fiuence of inhomogeneous couplings by a variety of ap- 
proximate numerical and analytical methods. Obviously 
the use of exact Bethe Ansatz methods to the general dy- 
namic problem would be a great advantage, but so far the 
possibilities have been limited to certain non-equilibrium 
situations in the related BCS-modeP^ and to the special 
case of a maximal bath polarization.^" Here we demon- 
strate how Bethe Ansatz results can be used to obtain the 
central features of the dynamic behavior for more general 
polarizations and coupling constants and compare with 
numerical complete diagonalization results. 

We first want to consider an initial state \L) = | JJ-i t? t 
, 4-, . . . , t)i where the central spin JJ- and Mb bath spins at 
specified sites L = {£i, . . . ,iMb} are in the down state. 
All other spins are in the up state. Hence \L) is an 
eigenstate of all Sf operators with total magnetization 
S^^^. = Nb/2 — All, — 1/2 and is initially not entangled in 



any way. The time evolution is given in terms of the 
eigenstates \M^) and eigenvalues Km^ of the model (1) 



\m) = Ee-'^-^*|Af.)(M.|L) 



(3) 



From (3) an explicit expression for the reduced density 
matrix of the central spin can be derived, ^"^ which we 
employ to evaluate 



{S^,){t) = \\l-2Y,\aj{tf 



with 



aj{t) = {J\Lit)) =J2{J\M.}{M,\L}e 



-itA, 



(4) 



(5) 



where both \L) and \J) are eigenstates of all Sf with the 
central spin fixed in the down state. Thus it suffices to 
specify only the Mb flipped bath spins in the subsector of 
dimension C^^ = iVf,!/((iV6 - Mb)lMbl). In Ref. [20] the 
matrix elements {M^\L) were explicitly given in terms 
of the quantum numbers (Bethe roots, see below) of the 
energy eigenstates \M,y) and of the 5f eigenstates \L); 
alternatively, they can be obtained from a complete di- 
agonalization. However, an exact calculation of the large 
number of terms in the sums (4) and (5) is impossible 
already for modest system sizes, except for special cases. 
In particular, a fully polarized bath S^^^ = Nb/2 — 1/2 
(i.e. Mb = 0) was studied in Refs. [10,20], where the sum 
(5) only contains Nb + 1 terms. 

We now would like to consider a more general polar- 
ization, Mb ^ 0, and single out the most important con- 
tributions to the sum (5). In this way, it is possible to 
estimate the dominant frequency scales and the widths 
of peaks in the frequency spectrum of Eq. (4) and thus 
to obtain the decoherence time. Our strategy is based on 
results in Ref. [21], where it was found that only a few 
product states \J) have an appreciable overlap (J|M^) 
for a given [M^) as long as Mb <C N. These product 
states are essentially those obtained from the classical 
ground state | |).,t:t7 •■•it) by flipping certain individual 
nuclear bath spins as outlined below, that spin pattern 
being also reflected in the local expectation values {Sj).^^ 
We illustrate this method for Mb = 1 first and generalize 
the results afterwards. 

The eigenstates \M^) can be classified by a set oiMb+l 
Bethe roots {wo,,^, . . . , cjMt,v} of the exact solution. Their 
positions in the complex plane are determined by coupled 
non- linear equations. ^^'^"^ The eigenvalues are given by 
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NbXi 



(6) 



fc=0 



To each eigenstate belongs a distinct root pattern that is 
related to the flipped spins relative to the all up state. ^^ 



In particular, a Bethe root in the origin corresponds ap- 
proximately to the application of a global lowering op- 
erator S't~t, a root Ag^i < lo < Ai induces essentially a 
superposition of states with spin flips on sites £, £ + 1^^ 
and a root ui = 0{Nb) mainly causes a flip of the cen- 
tral spin, respectively. Therefore, for Mb = 1 the state 
with the central spin and a single bath spin at site £ 
in the down state is most strongly overlapping with the 
six eigenstates that are characterized by two Bethe roots 
as follows: |0,0), |0,a;i_£_i), \0,oji,e), |0,wo), \uj[j,,uJo.i), 
\uj[ ^_^,ujQ£_i), where the roots can be approximately 
determined from the distribution of Aj in an expansion 
of rf := {NbXi)-^ and yi := dEf=i^? as^i 



Wo = l/d + yi + 0{d). 
ujo.e = 1/rf + yi - 2Ae + 0{d) 
uji^i « oj[g= A( + 0{d) 



(7) 



The sum of the squared overlaps |(Mi^|L)|^ from only 
those states yields 0.71 for £ = 1 increasing to 0.96 for 
£ = Nb^^ in a system with Nb = 15 and xi = i? = 2 in 
Eq. (2). Therefore, most of the weight in the expansion 
of \L) into eigenstates is indeed found by only consider- 
ing the six states listed above. For more homogeneous 
couplings (B=0.4) the corresponding overlaps are signif- 
icantly larger. 

Once these most important contributing states are 
known, the actual values of the overlaps {Mjj\L) are 
secondary, but the differences in the corresponding 
eigenvalues Am^ determine the spectral distribution in 
Eqs. (4) and (5) and therefore the decay time. Tak- 
ing the corresponding differences of eigenvalues using 
Eqs. (4-7), we find that the high-frequency contribu- 
tions occur at fio = ^{1 + dyi) and in an interval 
[J7i £ — A^/2, r^i £ -I- A^/2] around the ^-dependent fre- 
quency U.i^t = 2^ {l+dyi-2dAi), where A^ := (Af_i - 
Ai+i)/2 up to terms of order 0{d). The long time scale 
r^ 1/A^ resulting from the width of the peak near J7i ^ 
can accordingly be interpreted as the decoherence time 
for this dominant oscillation. In addition, there are low 
frequency contributions around ^2,f. — A^/2 -f 0{d), due 
to the inhomogencity in the couplings Aj^ which disap- 
pear for a homogeneous model. ^^ 

Before turning to other polarizations, Mb > 1, we 
can now also consider more realistic initial states with 
bath configurations other than the Sf product eigen- 
states hitherto considered. In the following we will use 
a uniform distribution over all bath states with a given 
polarization as the initial state, given by the density ma- 
trix 



A = ( cZ\ 



I^)(^|1a^„ 



(8) 



where Ijvfb is the projection operator onto the sector with 
Mb flipped bath spins relative to the fully polarized up 
state. In this case the frequency spectrum of (>5'o)(f) 
consists of the superposition of the spectra obtained for 
all allowed individual product states. Thus for Mb = 1, 
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FIG. 1: Time evolution (S'o)(i) (a,c) and the corresponding 
Fourier transform (b,d) for Nb = 15, xi = 2, and an initial 
uniform distribution of bath states with M(, = 1 (in a,b) and 
Mb = 2 (in c,d), respectively. The effect of two different 
homogeneity parameters B = 0.4 (circles, black) and B = 2 
(crosses, red) for the couplings in Eq. (2) is shown. Dashed 
lines (green) show the estimates for the decoherence times. 
The dotted lines (blue) give the position of the peaks in the 
homogeneous limit Aj =2. 



all frequencies ^ij, ^2,e contribute, with £ = 1, . . . , Nf,. 
Thus the high-frequency spectrum now consists of one 
sharp peak at Qq and of a second peak at f^i := fio ~ 2:1 
with a width A := Ai — Aj^^. Accordingly, the corre- 
sponding decoherence time for this oscillation is given by 



27r/A, 



(9) 



which is therefore directly linked to the inhomogeneity 
of couplings. Another way of interpreting these results 
is to say that the central spin precesscs in the effective 
field from the coupled bath state with a relaxation time 
Ti ^ T. The low-frequency tail contains frequencies in 
the interval [Aat,^, Ai]. 

In order to illustrate the decoherence process we show 
the time evolution of {S§){t) and the corresponding 
Fourier transform {Sq){uj) for Nt, = 15 bath spins as 
obtained from complete diagonalization in Fig. 1. Wc 
choose xi = 2 and two different values for B, correspond- 
ing to relatively homogeneous (B ~ 0.4) and relatively 
inhomogeneous (B ~ 2) couplings in Eq. (2). For nearly 
homogeneous couplings, the broadening of the peak near 
ill is demonstrated nicely. For B = 2, the two peaks can- 
not be distinguished any longer due to the large broad- 
ening of the peak around fli. 

We can generalize the above discussion to larger Mi, 
with the analogous initial density matrix in Eq. (8). For 
Mh = 2 three peaks are present centered around i7o and 
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FIG. 2: The time evolution (S'o)(i) (left) and its Fourier trans- 
form {S§){u}) for Stot = from complete diagonalization, for 
A'^i, = ll, 3^1=2, and two different choices of the homogene- 
ity parameter B in Eq. (2). The dotted lines (blue) give the 
position of the peaks in the homogeneous limit A = 2. 



fli as given above, and around fl2 '■= ^a — 2a;i. Neglect- 
ing complex string solutions of the Bethe Ansatz equa- 
tions and interactions between excitations, one expects 
the peak near i^2 to have twice the width A2 = 2 A. 
Again, in the homogeneous limit, the known^*^ frequen- 
cies A{Ni, -I- 1 — 2fc), k = 0,1,2, are recovered. From 
the beating of oscillations within the corresponding fre- 
quency ranges we estimate the overall decoherence time 
to be t' = 2r/3. The results for Mh = 2 in Fig. 1 are 
consistent with the estimates for r and t'. At the same 
time, one notices that additional spectral weight develops 
at small frequencies for increasing Mi, for inhomogeneous 
couplings. This can be traced back to the combination 
of states with Bethe roots (0,0;^^ ^) and (0, 0) in the sum 
in Eq. (4). Therefore, the smallest of these resulting 
frequencies is given by the most weakly coupled spins, 
fl ^ Ani^/2, around which indeed most of the spectral 
weight in the low-frequency region of Fig. 1-d develops. 

It is possible to increase Mi, further and continue this 
analysis with more peaks for other initial bath configura- 
tions as long as Mi, <C Ni,. The positions and widths of 
the peaks directly reflect the initially flipped bath spins 
and the choice of couplings. For larger Mi,, however, 
interactions between the elementary excitations may dis- 
tort the simple analogy between roots and flipped spins 
above. Therefore, we analyze the situation in the fol- 
lowing for small bath magnetization Mi, ^ Ni,/2 using 
complete diagonalization on a smaller system in order 
to see which features of the analytic considerations can 
survive in that case. 

In Fig. 2, we show results for A^f, = 11 in the subsec- 
tor S'fot = 0, where we took the average over all allowed 
product states as in Eq. (8) with the central spin point- 
ing down. For a weak inhomogeneity, one can still clearly 
distinguish the discrete peaks in the Fourier transform, so 



that the analytical predictions are still useful for M^ = 5 
in this case. The decoherence times of the lower fre- 
quency oscillations are again generally shorter, as can be 
seen by the widths of the corresponding peaks. 

This structure is lost for couplings with a significant 
degree of inliomogeneity, i.e. when the difference in cou- 
plings becomes larger than the peak separation. How- 
ever, as can be seen both in the time evolution and its 
Fourier transform, a new time scale at small energies oc- 
curs, which leads to low-frequency oscillations. In the 
Fourier transform, this shows up as a relatively strong 
peak at low frequency. It is reasonable to expect that 
this stems from the same low-frequency mechanism as 
discussed above, namely the overlap between the ground 
state - with all roots in the origin except for wo - and the 
lowest excited state - where one root is shifted out of the 
origin into the interval between AjVi, and Aat^-i. This 
yields dominant frequencies corresponding to the most 
weakly coupled spins, $7 ~ Ai^^/2^ leading to a charac- 
teristic long-time oscillation Att/An^ ■ The physical inter- 
pretation is that for generic disordered couplings fluctua- 
tions occur on all time scales up to the highest frequencies 
^0 ~ X^i ^j leading to a correspondingly large l/Ti, but 
for longer times there remains a relatively coherent low- 
frequency oscillation, which emphasizes the importance 
of the weakly coupled nuclear bath spins in the time evo- 
lution of the central spin.^^ This long time behavior can 
be explained by realizing that the most weakly coupled 
spins simply play the role of a relatively stable backgound 
field, but never actually become strongly entangled with 
the rest of the system. 

In summary, wc have analyzed the decoherence in the 



commonly used central spin model using exact methods. 
For relatively homogeneous couplings the locations and 
widths of dominant oscillations in the Fourier spectrum 
(5q)(w) can be predicted analytically even for smaller 
polarizations (large Mi,). The positions and widths of 
the peaks directly reflect the initial bath state and the 
choice of couplings. High frequency oscillations have the 
longest decoherence times. For larger inhomogcncity in 
the couplings the decoherence times become shorter and 
eventually the simple peak structure is lost. However, in 
this case the appearance of a low frequency feature can 
again be inferred from the root structure of the exact so- 
lution. Accordingly, it is possible to identify the physical 
behavior of the dynamics in the different cases, namely 
fast oscillations with decoherence due to the difference 
in coupling strengths for nearly homogeneous couplings 
on the one hand and relatively stable long time oscil- 
lations due to the most weakly coupled background for 
inliomogeneous couplings on the other hand. The results 
also provide a direct link from clear signatures in the dy- 
namics to individual excitations and characteristics of the 
model. This opens new possibilities and insights for the 
analysis of the dynamical behavior that can be obtained 
from independent methods, e.g. from more advanced nu- 
merical or even experimental studies. 
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